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CU ABSTRACT 

' The assumption of additive absorptivity by different components in compound particies 

^{■^ , is a wideiy used method applied in the literature to the analysis of the chemical 

£SJ ' and structural properties of astronomical (circumstellar, interstellar, protostellar, and 

cometary) silicates as well as other materials. The errors intrinsic in this additivity 
assumption, which, in application to astronomical spectra, amounts in some case to 
spectrum subtraction have not always been adequately considered in previous works 
' on silicate mineralogy. The failings in the "spectrum subtraction method" (intrinsically 

(f) , the same as the additive absorptivity assumption) are discussed here in terms of 

silicate core-ice mantle grains with various shapes. It is shown that these assumptions 
result in substantial errors for spherical grains. For spheroidal grains, the errors are less 
significant and the spectrum subtraction method can be used to remove the ice mantle 
effects. It is demonstrated that there is no significant improvement by considering a 
distribution of spheroidal shapes. It is further shown that the presence of additional 
organic mantles substantially modifies the silicate mineralogy interpretation. 

2 ' Key words: ISM: dust, extinction - infrared: ISM: lines and bands - scattering 

CO 

C3 

> '■ 1 INTRODUCTION 

• i-H , 

' Silicate grains are ubiquitously present in various cosmic environments. Theoretical efforts have been continuously made to 
^ , identify the detailed chemical and structural properties of astronomical silicates (Draine & Lee 1984; Ossenkopf, Henning, & 
Mathis 1992; Greenberg & Li 1996; Mathis 1998; Demyk et al. 1999). On the other hand, during the last decade there has 
been considerable progress in laboratory measurements of astronomical silicate analogues (see Dorschner 1999 for a review). 
More recently, the high quality spectroscopic observations in the 10, 18/im (as well as the longer wavelength range) silicate 
bands region of a wide range astronomical objects by the Infrared Space Observatory (ISO) together with the rich laboratory 
data set offer us an unprecedented opportunity to study silicate dust mineralogy. 

The analysis of astronomical silicate composition is commonly based on the direct comparison of astronomical absorp- 
tion/emission features with laboratory data of terrestrial silicate candidates or with theoretical spectra calculated from the 
optical constants measured for silicate analogues. However, in most regions silicates are combined with other materials. For 
example, Greenberg (1978) proposed that in the diffuse interstellar medium silicates are coated by a layer of carbonaceous 
organic refractory; whereas Mathis & Whiffen (1989) suggest that interstellar silicates are combined together with amorphous 
carbon grains to form porous composite aggregates. GEMS (silicate glass with embedded metal and sulfides) in interplanetary 
dust particles (IDP) which provide a spectral match to the interstellar 10/im Si-O stretching band is usually mantled with 
or embedded in amorphous carbonaceous material (Bradley et al. 1999). Cometary dust particles are also believed to be a 
mixture of various forms of silicate and organic ("CHON" particles) materials (Greenberg 1982; Kissel et al. 1986; Wooden 
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et al. 1999). In dense molecular clouds or protostellar disks, the consensus view is that refractory grains are further coated by 
an outer ice mantle. 

In studying silicate composition, it is thus critical to accurately remove the effects caused by the grain's heterogeneous 
nature. A simple and commonly adopted method is to simply subtract the non-silicate contributions (determined from their 
characteristic absorption/emission properties) from the entire observational spectrum and then make a comparison between 
the residual spectrum either directly with laboratory data or with a theoretically derived absorption spectrum by particles 
using the laboratory optical constants. Hereafter we term this the "spectrum subtraction" method. 

For example, in a detailed effort to infer the precise composition of the silicate dust around two massive protostellar 
objects RAFGL 7009S and IRAS 19110+1045, Demyk et al. (1999) subtracted the ice absorption from the total observed 
absorption to obtain what they presumed to be the observed silicate absorption spectrum as if it could be obtained by adding 
"ice" absorptivity to the absorption spectrum for bare spherical silicate material grains. Two key assumptions here, both of 
which affect the derived absorption, are that the grains are spherical and that the absorption by a core-mantle particle is the 
same as the absorption by the core plus the absorption by the mantle. We term this the "absorption additivity assumption". 
Note the spectrum subtraction method and the absorption additivity assumption are intrinsically the same. On the basis of 
this method, Demyk et al. (1999) then modelled the observed silicate absorption spectrum by performing Mie calculations 
(together with effective medium theory) of spherical porous grains of various compositions. 

In the current studies of the crystallinity of circumstellar silicate grains, the circumstellar dust emission spectra are often 
modelled as co-added emissivities of individual dust components in spite of the fact that (at least) some dust components 
are actually not separated (e.g., Malfait et al. 1998). The absorption additivity approach has also been applied to cometary 
silicate analysis (see e.g. Brucato et al. 1999). 

The absorption additivity (spectrum subtraction) method takes as a premise that the coupling between the silicate and 
non-silicate materials is negligible. It is obvious that such a premise is physically incorrect (e.g. see Bohren & Huffman 1983). 
However, this fundamental question has not always been adequately considered in modelling the silicate emission/absorption 
features in the literature. An equally important problem has to do with the effect of particle shape on the absorption spectrum. 

The purpose of this work is not to model any specific astronomical objects, but to address the question of to what extent 
or in what circumstances such a premise is valid or acceptable. In doing so we will consider a simple case: silicate grains coated 
by an ice mantle as frequently proposed in the literature although it is likely that an inner carbonaceous mantle may also be 
present (Greenberg 1978). In §2 we summarize the optical constants employed in this paper. In §3 we discuss the spectrum 
subtraction method in terms of spherical shape. §4 is devoted to spheroidal shape. We extend then to a distribution of shapes 
in §5. In §6 we discuss the implications of a more realistic model; i.e., the coated core-mantle spheroidal grains. We summarize 
our results in §7. In a subsequent paper we will perform detailed modelling on the interstellar silicate bands and discuss the 
depletion of heavy elements and its implication on dust models (Zhao, Li, & Greenberg 2001). 



2 OPTICAL PROPERTIES OF DUST MATERIALS 

We consider three types of dust materials: amorphous silicates, carbonaceous organic refractories, and ices. We take the 
dielectric functions of amorphous olivine MgFeSi04 and organic refractory materials from Li & Greenberg (1997). For ices, 
we consider a mixture of H2O : CH3OH : CO2 = 1.3 : 1 : 1 (80-90K) which provides the best interpretation of the ISO spectra 
of RAFGL 7009S (Dartois et al. 1999). The incorporation of CH3OH and CO2 adds substantial ice mantle effects to the 
silicate features since they have strong absorptions such as the 9.75^im C-O stretching mode of CH3OH, and the 15.2/im C-O 
bending mode of CO2. We note that neither the specific choice of silicate optical constants nor the precise constituents of the 
ice mixtures will affect our conclusion since in the present paper, we are not intending to reproduce any specific observations 
but rather to demonstrate the effects of grain shape and the presence of an ice mantle and therefore to explore the validity 
of the spectrum subtraction method. We defer our detailed modelling of astronomical silicates to Zhao et al. (2001) where 
the chemical variations of silicate, carbonaceous, and ice materials, and grain shapes and morphology (porosities) will be 
discussed in detail. 

Since there are no existing optical constants available for such an H2O : CH3OH : CO2 = 1.3 : 1 : 1 mixture, we construct a 
set of synthetic refractive indices from the existing data for pure H2O at 80K, pure CH3OH at 75K, pure CO2 at 70K (Hudgins 
et al. 1993; with mass densities » 1.00, 1.23, 1.56 gem -3 , respectively) by applying the Bruggeman Effective Medium Theory 
(EMT) (Bohren & Huffman 1983) and assuming a volume fraction of H 2 : CH3OH : CO2 = 1:1.1: 1.2. Theoretical spectra 
calculated from the synthetic optical constants are consistent with the laboratory IR spectrum of Ehrenfreund et al. (1999). 
Again, we emphasize here that the precise knowledge of the optical constants of ices is not critical since the purpose of this 
paper is not to infer the exact dust composition. 
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3 SPHERICAL GRAINS 

We first consider spherical grains. Under the condition of x — 2na/\ <C 1 and \mx\ <C 1 (a is the grain radius, m is the 
grain refractive index which relates to the dielectric function e through e = m 2 ), the absorption cross section of a core-mantle 
spherical grain C a bs can be calculated from the Rayleigh approximation (van de Hulst 1957) 

c»./v = 3 ^im ( h - ;;s e - ± ± if:- ± ^ - 1 ( d 



A [ (e m + 2)(e c + 2e m ) + /c(2e m - 2)(e c - e m ) 

where V is the dust volume; e c and e m are the complex dielectric functions for the core and mantle materials, respectively; f c 
is the volume fraction of the core. The scattering cross sections C sca are much smaller than Cabs so that the extinction cross 
sections C ex t ~ C a bs- 

We carried out calculations for spherical silicate core-ice mantle grains, taking the silicate core size a s u = 0.07/im, a 
typical size for the interstellar silicate dust (Li & Greenberg 1997). Note the knowledge of the exact grain size is not critical 
since interstellar grains (~ 0.1/im) are in the Rayleigh limit in the wavelength range of interest here and also because of this 
no size distribution is needed. 

In Fig. la we plot the absorption cross section of the core-mantle grain 

C tbt iCe with fflsii=0.07^m, / si i=2/3 (the volume 
fraction of the silicate core, si — Vicc/Usii+ice) in the silicate features region. Also plotted are the absorption cross sections 
for a pure silicate grain C a b a with the same size as the silicate core, and from a pure ice grain C a b C s with size ai cc ~0.0556/im 
corresponding to the volume of the ice mantle. Note that the pure ice grain spectrum closely reproduces the ice features 
of the core-mantle grain. According to the spectrum subtraction method, the contribution of the silicate core to the total 
absorption cross section « C° b '~ SSM = C'° b 1 s '" lcc — C a b C s - However, as clearly illustrated in Fig. la, the C^ b '~ SSM spectrum deviates 
significantly from that of the pure silicate grain C a b s : the peak position of both the 10pm Si-O and the 18/im O-Si-O features 
shift to longer wavelengths; both features have a broader red wing with an error xssm = /C^tm — 1 ~ 8% for the 

10/im wing and ~ 12% for the 18/^m wing. These differences are comparable to some of the differences ascribed by Demyk et 
al. (1999) to different materials and/or different porosities (for example, the simple spherical enstatite and olivine fits to the 
silicate absorption spectra [with ice absorptions subtracted] of RAFGL 7009S and IRAS 19110 which were rejected by Demyk 
et al. fail to reproduce the red wing by < 10% for A < 18/^m and < 20% for 18 < A < 20zjm). In light of this, Demyk et 
al. (1999) 's invoking of an admixture of iron, aluminum and some degree of porosity to explain the red wing of the 18/^m 
feature seems unnecessary since it may be just an artifact of the spectrum subtraction method they adopted (see their §4.2 
and Fig. 4). 

We have also considered the variation of the ice mantle thickness. For illustration we present in Fig. lb the errors xssm 
caused by the spectrum subtraction method for /i cc = Vice/Vsii+ice =1/5 (dotted), 1/3 (solid), and 1/2 (dashed), respectively. 
Fig. lb shows that the spectrum subtraction method can not give acceptable results even for a very thin layer of ice mantle 
(f«0.0054 M m for / ice =1/5). 



4 SPHEROIDAL GRAINS 



Cosmic dust grains are unlikely to be spherical. On the contrary, polarization of starlight implies both nonsphericity and 
alignment. We now consider spheroidal shape which can be exactly solved in the Rayleigh limit. For an ensemble of randomly 
oriented spheroidal grains^, the absorption cross section is (Greenberg 1972; Lee & Draine 1985) 



C a bs — — C 



a bs ' g^abs 



(2) 



where C^ bs and C^s are the absorption cross sections for light polarized parallel and perpendicular, respectively, to the grain 
symmetry axis. In the Rayleigh limit, for confocal core-mantle spheroidal grains we can calculate the absorption cross sections 
Cabs from ( Gilra 1972 ; Draine & Lee 1984) 



Cabs/^flm. 



(em - 1) 


Lc ( e c ^m) ~\~ 


+ /o(ec - e m ) 


L|k X (l-e m )+e m 




Lc (Cc ^m) H~ £m 




i + 4^(^-1) 


+ / c L|1; X (1 - Ljk ± )(e c - e m )(e m - 1) 



(3) 



t Similar results are expected for aligned grains; for example, the deviation \' = y S { [^abs"' 06 — Tabs] /Tabs — ^} ^ or 
wavelength range of 8.5 < A < 20^tm (N is the number of data points) from 3:1 prolate grains of perfect spinning alignment with 
/icc=l/3 is about 0.023 in comparison with \' = 0.020 for randomly oriented prolates. 
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Figure 1. a: Absorption cross sections of spherical grains. Dashed (C*j J 1 + ioe ): silicate core-ice mantle grains with a silicate core radius 
0.07^tm and a silicate volume fraction / s ;i = 2/3 (i.e., /i cc = 1/3); long dashed the volume-equivalent sphere for the ice mantle; 

,sl1 ,): a pure silicate sphere with the same size as the silicate core of the core-mantle dust; solid: (^ s jl~ SSM j n the framework 



of the SSM, one should have C 



il-SSM 



abs ' 



The contrast between C'\ 



il-SSM 



dl-SSM 
abs 

, ~ ^ a b s - >■"*> wutiaoi/ ucuwra] ^abs (solid) with CfjL (dotted) indicates the unapplicability 

of the SSM to spherical grains, b: The errors XSSM caused by the spectrum subtraction method for /; cc =1/5 (dotted), 1/3 (solid), and 
1/2 (dashed). The errors in the red wings of the 10/xm and 18^m features are comparable to the differences ascribed by Demyk et al. 
(1999) to different materials and/or different porosities (see our text and their §4.2 and Fig. 4). The sharp features at 3.1, 4.2, 9.75, and 
15.2/^m are due to the incompleteness in subtracting the ice features. This will not affect our conclusion since in this work we are mainly 
interested in the broad 10 and 18^tm silicate bands. 



where L\, L c , L'm, and L m are the depolarization factors of the core and mantle parallel and perpendicular, respectively, 
to the grain symmetry axis. The depolarization factors i", L x relate to each other through L x = (1 — -£/")/2 because of 
rotational symmetry and characterize grain shape via Eq.(4): 



^£ [^ ln (f^f ) - !] . for P rolate ( a > & ) 

(l - Itan^e) , for oblate (a < b) 



(4) 



where the eccentricity e is given by the two semiaxes a and b (with a semiaxis along the symmetry axis and b semiaxis 
perpendicular to the symmetry axis) 
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_ J y/l — (b/a) 2 , for prolate (a > 6) 
6 ~ \ ^(b/a) 2 - 1, for oblate (a < 6) 

and the eccentricities of core and mantle e c , e m (e c > e m ) relate to each other through 
f ~S S ~ T~T2~ i f o r prolate (a m > 6 m ) 

/« = j: . M w (6) 

[ "1^I+S"> for oblate (a m < b m ) 
In the special case of e c = e m = e, Eq.(3) reduces to the solution for a homogeneous spheroidal grain 



A \(e-l)LW^ + 1 

By applying the method described above, we calculate the absorption cross sections of confocal silicate core-ice mantle 
spheroidal grains with a volume-equivalent sphere radius (the radius of the equivalent sphere having the same volume as the 
spheroid) r^J = (ab 2 ) 1 ^ 3 = 0.07/im for the silicate core and a variable ice mantle controlled by / s n. We first consider a/b = 3 
prolate core-mantle grains with r2j = 0.07/im (a/b = 3 is for the ice mantle; the corresponding a/b for the silicate core is 
ft! 3.57 if / s ii=2/3). The reason for choosing an a/b = 3 prolate shape is that Greenberg & Li (1996) found that it provides the 
best fit to the 10, 18/im silicate polarization features of the Becklin-Neugebauer (BN) object in terms of ice-coated core-mantle 
dust grains. Also calculated are 1) C'^ s , the absorption cross section of the volume-equivalent pure ice prolate (a/b = 3); 2) 
CJbsi the absorption cross section of the volume- equivalent pure silicate prolate with a/b » 3.57 determined from Eqs.(5,6); 
and 3) the ice mantle "subtracted" silicate absorption C°5J7 SSM - in Fig. 2a we plot the corresponding errors xssm for /i Ce =1/5 
(dotted), 1/3 (solid), and 1/2 (dashed). 

As can be seen in Fig. 2a, C^ 1 ~ SSM is closer to C^ B in comparison with the spherical model (§3). The errors xssm for 
the red wings of the 10/im and 18/im features are respectively ~ 2.5% and ~ 2% (for / s n=2/3), in sharp contrast to those of 
the spherical model (~ 8% and ~ 12%; see Fig. lb). 

Similarly, we carried out the same calculations but for oblate grains with a/b = 1/2 (the core has a/b « 1/2.52). The 
a/b =1/2 oblate shape was chosen since it was shown to provide a better match than any other shapes to the 3.1/im ice 
polarization feature (Lee & Draine 1985) and the 10pm silicate polarization feature (Hildebrand & Dragovan 1995) in terms 
of ice-coated silicate grains (without an intermediate carbonaceous mantle). In Fig. 2b we plot xssm for /i cc =1/5 (dotted), 
1/3 (solid), and 1/2 (dashed) for oblate grains. The results are very similar to the a/b = 3 prolate model (Fig. 2a). The errors 
are < 3% for the red wings of the 10/mi and 18/im features. 

From these results we conclude that for spheroidal grains, the spectrum subtraction method works better than for spherical 
grains^]. Lacking prior knowledge of the true shape of dust grains, the spectrum subtraction method, as a first approximation, 
can be adopted to study the silicate mineralogy if grains are taken to be spheroidal. This does not necessarily imply that 
grains need to be spheroidal. What it does imply is that the errors caused by the spectrum subtraction method, in reducing 
the silicate mineralogy in terms of a spherical shape, are larger than those of spheroidal shapes. If the choice of spheroidal 
shapes seems rather arbitrary despite the fact that the error xssm is indeed smaller, we remind the reader that the spectrum 
subtraction method is physically incorrect. A better approach to the silicate mineralogy problem is to carry out iterative 
modelling of the observational spectra using Mie theory together with the effective medium theory or the discrete-dipole 
approximation for various shapes (e.g. see Zhao et al. 2001). 



5 GRAINS WITH A CONTINUOUS DISTRIBUTION OF ELLIPSOIDS (CDE) 

The shape of a spheroidal grain is characterized by its eccentricity e (and the depolarization factors L", L x via Eq.[4]). We 
now extend §4 where only a single shape (a single eccentricity) is considered to an average of a distribution of shapes (for the 
same volume grain). We assume two kinds of shape distribution functions: 1) dP/dL^ = constant, i.e., all shapes are equally 
probable (Bohren & Huffman 1983); 2) dP/dL n = 12L"[1 - L 11 ] 2 (Ossenkopf, Henning, & Mathis 1992; Will & Aannestad 
1999) which peaks at spheres (L" = L x = l/3) and is symmetric about spheres with respect to eccentricity and drops to zero 
for the extreme cases (e — > 1 [L" — > 0] infinitely thin needles or e — > oo [L" — > 1] infinitely flattened pancakes). Averaging 
over the shape distribution, we have the resultant absorption cross section Cabs = J Q dI/"dP/dL"C a bs(i") where C a bs(i") is 
the absorption cross section of a particular shape L" (note L" is for the mantle; the core depolarization factor is derived from 



t Why the spheroid model works better is not obvious but it may be related to the fact that the spheroidal absorptivity follows the 
"material" absorptivity (m ) while spherical one does not (Greenberg et al. 1983; van de Bult et al. 1985). 
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Figure 2. The errors XSSM caused by the spectrum subtraction method for /i C e=l/5 (dotted), 1/3 (solid), and 1/2 (dashed) - a: for 
a/b = 3 prolate; b: for a/b =1/2 oblate. In comparison with the spherical basis, the spectrum subtraction method based on spheroidal 
grains provides a better approach. The errors at the red wings of the 10^m and 18/im features are much smaller than those of the 
spherical model (see Fig. lb). 

Eqs.[4-6]). Again, we assume confocal geometry for core-mantle grains, with the above dP/dL^ (as well as e) applying to the 
outer surface. 

We calculated the averaged absorption cross sections of grains with a continuous distribution of ellipsoids (CDE) with 
r^=0.07/xm, / s ii=2/3 for dP/dL ]l = constant and dP/dL 11 = 12L 11 [1 - L 11 ] 2 respectively. We also calculated (1) C^ s , the 
averaged absorption cross sections predicted from the silicate core itself with eccentricities obtained from Eq.(6) for a given 
set of mantle eccentricities e and f B n; (2) C^bs, the averaged absorption cross sections of ice spheroids with the same volumes 
and the same eccentricities e as the ice mantles of the core-mantle spheroids; and (3) C'^J SSM . 

It is found that the C°j 3 '~ SSM of the dP/dL^ = constant model gives similar results to those of the single-shaped spheroidal 
grains (§4) and leads to smaller errors than the spherical model (see Fig. 3a). In contrast, although the dP/dLfi = 12L" [1 — L"] 2 
model was generally considered as "more physically reasonable" (e.g. see Ossenkopf et al. 1992), it results in a much larger 
error (even for /i ce = 1/5, see Fig. 3b). In particular, C^" 881 ^ deviates from CJ^s significantly in the peak regions of the 
10, 18/xm. This can be understood by the fact that the dP/dL^ = 12L"[1 — L"] 2 distribution peaks at spheres for which the 
spectrum subtraction method results in a significant error (§3). In comparison with §4, we conclude that the CDE model 
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Figure 3. The errors XSSM caused by the spectrum subtraction method for /i cc =l/5 (dotted), 1/3 (solid), and 1/2 (dashed) - a: 
for a continuous distribution of spheroids with all shapes equally probable (i.e. dP/dL^ =constant); b: for a continuous distribution of 
spheroids with the shape probability distribution dP/dL^ = 12lJ [1 — L^] 2 . The silicate core size is always taken to be r^. c }=0.07/im. 
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does not seem to be better suitable to be applied within the framework of the spectrum subtraction method than the single 
(non-spherical) shape model. 



6 COATED CORE-MANTLE GRAINS 

Finally, we consider a model in which the silicate core is coated by an carbonaceous mantle (in addition to its outer ice mantle) 
(Greenberg 1978). We emphasize here the fact that with or without ice mantles there is generally no pure silicate core but 
rather one which bears an organic mantle as seen in the diffuse interstellar medium, cometary dust, and as implied in molecular 
clouds (see Greenberg & Li 1999 for a summary). The absorption cross sections for an equal- eccentricity three-layered silicate 
core-carbonaceous inner mantle-ice outer mantle grain is (Farafonov 2001; Voshchinnikov & Mathis 1999) 

27t47T t f A}^ - A\ 

3 1(4'" 

© RAS, MNRAS 000, [!]-?? 



8 A. Li, J.M. Greenberg, & G. Zhao 



( A- 1 - \ / l ifl> x \ ( (c ijn /c ro -i)L |l - L + i (eWe™ - l)^ 1 '^ 11 ^ - l)//™ A 

^ 4' X / " V e ° m e om (i ll ' ± -l) y ^ -(e im /e om - 1)/™ -(«W««» - ±){L h± - 1) + 1 ) , , 

f (ec/e im - 1)L^ + 1 \ W 

^ -(Cc/Cim - l)/e j 

where e c , fim, £om are the dielectric functions for the silicate core, the organic refractory inner mantle, the ice outer mantle, 
respectively; f c and /i m are the volume fractions of the silicate core and the organic refractory inner mantle, respectively. 
Note that, in contrast to §4, the depolarization factors here are the same for all layers. 

We have calculated the absorption cross sections of ice-coated silicate core-organic refractory mantle grains with r° ;l =0.07pm, 
/sii = /or=2/5, /ice=l/5, for a/b = 3 prolate as well as for a/b = 1/2 oblate, where f s u (= / c ), f OI (= ,fi m ), and /i ce 
(= 1 — /o — /im) are respectively the fractional volumes of the silicate core, the inner organic refractory mantle, and the 
outer ice mantle. As not unexpected, C°J 3 1 1 7 SSM (= C'° 1 b 1 + or+lcc — Cabs) is very different from C^L in the peak positions, the 
feature widths, and the feature strengths of both the 10/im and 18/im features. Since the organic refractory material has little 
absorption around that region, the inclusion of such an inner mantle significantly suppresses the 10/j.m feature (not shown 
here). It is obvious that silicate grains with carbonaceous coating will not behave like naked silicates. Therefore, it is essential 
to take into account the possible heterogeneous nature of cosmic dust when modelling the silicate mineralogy. 

7 SUMMARY 

We have investigated the validity of the "spectrum subtraction" method or the "absorption additivity assumption" widely 
adopted in studies of astronomical silicate mineralogy when interpreted using silicate core-ice mantle grains with spherical, 
spheroidal, and a continuous distribution of ellipsoids (CDE) shapes. We found that: 1) the spectrum subtraction method 
results in significant errors for spherical grains; 2) for spheroidal grains, the errors are much less significant and the spectrum 
subtraction method is reasonably successful in removing the ice mantle effects; 3) grains with a uniform distribution of 
spheroidal shapes produce results little different from those for a single eccentricity; 4) as expected, grains with sphere-peaked 
distribution of spheroidal shapes, lead to appreciable errors. In addition, we have also discussed the effects caused by the 
inclusion of an intermediate carbonaceous mantle which is often not taken into account in silicate considerations. As not 
unexpected, the spectrum subtraction method applied here leads to even more unacceptable errors. 

Our conclusion is that, in order to correctly infer the precise composition of astronomical silicates, care should be taken 
to appropriately take into account the heterogeneous (chemical and/or structural) nature of dust grains. The spectrum 
subtraction method should be used together with the non-spherical (e.g., prolate or oblate) assumption. It does not appear 
to have obvious advantages to use a distribution of spheroidal shapes. 
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